Multi-level, multi-party singlets as ground states 
and their role in entanglement distribution 
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We show that a singlet of many multi-level quantum systems arises naturally as the ground state 
of a physically-motivated Hamiltonian. The Hamiltonian simply exchanges the states of nearest- 
neighbours in some network of qudits (d-level systems) ; the results are independent of the strength of 
the couplings or the network's topology. We show that local measurements on some of these qudits 
project the unmeasured qudits onto a smaller singlet, regardless of the choice of measurement basis 
at each measurement. It follows that the entanglement is highly persistent, and that through local 
measurements, a large amount of entanglement may be established between spatially-separated 
parties for subsequent use in distributed quantum computation. 
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Entanglement between spatially-separated systems is a 
pivotal resource in quantum information theory, enabling 
distributed or networked quantum computation. Thus 
there is an enormous interest in extracting this^ resource 
from the ground states of many-body systems [1], in par- 
ticular through measurements 2]. In such schemes, the 
measurement bases have to be carefully optimised. Could 
there be other systems whose ground states offer a more 
flexible method of entanglement extraction, and could 
the amount of entanglement exceed the currently-known 
limits? Here we show that this is indeed the case for a 
particular system: a frustration-free permutation Hamil- 
tonian obtainable from the many-species Hub- 
bard model (either fermionic or bosonic) . Such Hamilto- 
nians have provoked much interest in many-body theory, 
since permutation Hamiltonians are a generalisation of 
the Heisenberg model [sj, and the Hubbard model may 
be used to model phenomena ranging from superconduc- 
tivity to ultracold atoms in optical lattices. We show that 
the ground state of this Hamiltonian is the iV-partitc sin- 
glet made of A^-level systems, independently of both the 
spatial arrangement of the subsystems, and the strength 
of their couplings. 

The singlet is a highly-entangled, curious state, and 
has been shown to have wide-ranging applications in 
quantum information: in problems lacking classical solu- 
tions , for multi-party remote state preparation j^l, and 
for encoding qubits in decoherence-free subspaces [3, . 
In this Letter, we further consider some of the state's 
properties in the context of quantum information. We 
show that if N parties share an A^-level singlet, and M 
parties perform successive measurements (each in a ran- 
dom basis) , the remaining parties share a singlet of N—M 
systems, regardless of the choice of measurement bases 
(measurement in the same basis has been previously con- 
sidered 0])- A direct consequence of this is that these 
states have the highest possible persistency of entangle- 
ment (the robustness of a multi-party entangled state 
to local measurements [icj ) . This basis-independence is 



particularly interesting in the context of entanglement 
distribution: as already stated, entanglement can be es- 
tablished between spatially-separated parties (for exam- 
ple, between the extremal spins of a ID spin chain) by 
performing local measurements on the other, interven- 
ing, systems. One can define the localisable entangle- 
ment as the maximum entanglement establishable in this 
way In general, maximising this quantity requires 
careful optimisation of the local measurement bases, and 
averaging over all outcomes. However, for the singlet, 
even if one selected the basis randomly at each measure- 
ment, one can always establish a maximally- entangled 
state between the unmeasured parties. This is a quali- 
tative difference in mechanism of entanglement localisa- 
tion to those systems studied previously 2\. Moreover, 
another important difference emerges if one considers the 
entanglement establishable between two subsystems {i.e. 
groups of constituent systems), rather than two individ- 
ual constituents, of a system {e.g. blocks of spins in a ID 
chain, rather than individual spins). In models studied so 
far, the entanglement establishable between two subsys- 
tems has no reason to differ from that of two individual 
parts; however, in our model two n-particle subsystems 
may share a maximally-entangled state of dimension (^J') 
(equivalent to log2 (^^) EPR singlets). We also show that 
the block entropy of an A^-singlet (the entanglement of a 
subsystem and the remainder of the system) is log2 (^) 
where L is the block size. Finally, we discuss a potential 
physical realisation of a singlet in an optical lattice. 

Qudit singlets: A qudit is a generic d-level system, and 
qudit singlets |<S'^'*) are A^-partite states with the prop- 
erty U'^^\S^^'^) = 15^^) (up to a global phase), where U 
is an arbitrary one-qudit unitary operation. For the case 
d — N , they may be written 
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is the generalised Levi-Civita symbol. 
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and the state is written in the basis {lai)}^^^ each 
qudit. The sum is taken over all permutations of in- 
cides {rii, • • • , npf}, and crucially, is anti-symmetric with 
respect to the exchange (permutation) of any two qu- 
dits: I.e. P,,\s'-pia)) = -\S'j^\a)) for i,j e [1,7V] 
where for given states |V')i G Hi,\(f>)j e Hj, the per- 
mutation operator swaps the states according to 
Pij\ip)i\4>)j = \4>)i\''P)j and has eigenvalues ±1. 

Permutation Hamiltonian ground states: The famil- 
iar isotropic Heisenberg Hamiltonian (cr* • cr^) acting on 
two qubits i, j is equivalent (when the two-qubit iden- 
tity operator is added) to the permutation operator . 
The natural generalisation of this Hamiltonian for d-level 
systems is a sum of permutation operators Py For 
a general network of qudits (some arrangement of qu- 
dits connected by permutation operators) one may use 
the language of graph theory and ascribe a finite graph 
G := {V{G),E{G)}, where F(G') denotes its set of ver- 
tices and E{G) its set of edges — if are adjacent ver- 
tices, G E{G). At each vertex we associate a d-level 
Hilbert space. We avoid the term lattice at this stage, 
since this implies a regular arrangement of qudits; these 
results hold for more general networks, making graph 
theory the natural description. The Hamiltonian is then 
written H — J^i jeE{G) '^u^y ' ^^^^ ^lay be obtained from 
the many-species Hubbard model (discussed below). We 
now show that for these Hamiltonians, with d — N and 
all couplings Jij > (antiferroniagnetic) , the qudit sin- 
glet |S'^'') arises as the ground state. We term this an 
W-singlet' [24]. 

Lemma 1. The lowest possible energy state of a permu- 
tation Hamiltonian has energy equal to that of an eigen- 
state of all {Pjj]i, j S E{G)} with eigenvalue —1. 

Proof. By definition, the ground state must minimise 
the energy (■0 1 1 ^/;) . Now minj^)g(cd)»jv > 
SijeB(G) '^i'^IV'>G(C'')®" (V'l'^y^'iilV') ■ The smallest 
eigenvalue of P^ is —1, thus min (V'iPy IV") — — 1; 
it follows that min|^)g(cd)®N (t^ ] -ff ] V') > -J2i,j'^ij- 
Equality exists for an eigenstate of all terms in the 
Hamiltonian {i.e. {Pij\i,j E E{G)}) and if this state 
exists, it is the ground state. □ 

Lemma 2. // a state is an eigenstate of {Pij\i,j G 
E{G)}, it is an eigenstate of all {Pij\i,j G [1, A^]}. 

Proof. Consider first a linear chain with nearest- 
neighbour permutations. Any other permutation may 
be written as a product of an odd number of nearest- 
neighbour permutations, {e.g. P13 — P23P12-P23; gen- 
erally Pi,fc = Pi,fc-iPA;-i,fc-Pi,fe-i). Thus an eigenstate 
of all nearest-neighbour permutations must also be an 
eigenstate of all possible permutations. This can be read- 
ily generalised to any connected graph, since any 'path' 
through the graph is equivalent to a ID chain of nearest- 
neigbour connections; thus any permutation in any graph 



can be written as a product of an odd number of nearest- 
neighbour permutations. □ 

Theorem 3. The ground state of a permutation Hamil- 
tonian on N N -level systems is an N -singlet. 

Proof. In the above lemmata we have shown that a state 
antisymmetric under all permutations is a valid ground 
state; since a qudit singlet for d = N hy definition satis- 
fies this, it is a valid ground state. It can easily be shown 
to be the unique ground state by assuming the existence 
of another distinct ground state \(j)) ^ ]S'^''). Being a 
ground state, \<j)) must minimise the total energy, leading 
to {(t)\Pij\<f>) = -1 for ah i,j. This implies ]</)) = JS"^'). 
This is a contradiction and thus there cannot be another 
ground state. □ 

Having shown how such singlets arise as ground states, 
we now go on to consider some of their properties and 
applications. 

Local measurements: We show that when A'^ parties 
share such a state, and some of them perform measure- 
ments on their qudits in randomly chosen bases (vary- 
ing from party to party), a qudit singlet is still es- 
tablished between those unmeasured qudits, in a ba- 
sis related to that used for the last measurement. To 
demonstrate this, we will make use of the property Q 
U ®t®^-^\S^^^) = I® C/t®^-i]5(f)) and introduce 
the notation \S^j^Si\(^] Pi)) to denote an (A^— l)-singlet 
written in the basis {l/3i)}^i at each qudit, with level 
]/?;) absent {i.e. Trpi]/9/)(/3;] = for all i where pi is the 
reduced density matrix of the zth qudit). 

Theorem 4. When an N -singlet \S'^^\a)) written in 
a basis {\ai)} is measured at one site using an arbitrary 
basis {\(3i)} :— {U\ai)}, the state obtained is a product 
of a state \Pi) at the measured site and a smaller singlet 
\sl^_i^\f3; Pi)) written in the basis {\Pi)}. 

Proof. Consider the outcome when we perform a von 
Neumann measurement {\Pi) at one qudit. With- 
out loss of generality, we take the above measurement on 
the qudit labelled 1. Since \(3i) = U\ai) (a^] we 
may write 

= C/]a.)(a.]C/t®I«^-V^^'(a)) 
-{/]a.)(a,]I®C/®^-Vf^(")>- (1) 

To proceed, we write the A^-singlet as 
1 ^ 

]4^)(a)) = ^^(-)^+V,)i]5f_",i)(«;a,))2.....Ar, 



3 



and it thus follows that 



i®c/®^-vr("))i,-,A^ 

N 



Whilst \s'}^Si'\oi; ai)) is a singlet written in the 
{|Q;;)}/^i basis at each qudit, one can see that the 
term U®'^~^\si^-^^\a; at)) is a singlet written in the 



basis [25|. Thus if the measurement outcome is 
the overall state is projected to 

-(-r+'IA)i®l^fri'^(/3;A)>2,...,^. 

This completes the proof. □ 

The significance of the above theorem is revealed when 
one considers successive measurements at different qu- 
dits. Indeed, by iterating the above proof, it becomes ap- 
parent that by measuring in a different basis at each qu- 
dit, the remaining, unmeasured qudits will be projected 
to a singlet in a basis related to that used for the final 
measurement. 

Corollary 5. If M parties perform successive measure- 
ments in arbitrary bases (where the mth party uses the 

basis Bra = {|af^)}.^i,..,™-i = {Y{T=iU^'\<^f)}), the 
remaining parties share an (N — M) -singlet in the Bj^i 
basis at each qudit, with the restriction that at each mea- 
surement, the basis transformation operates on a space 
whose dimension is one less than the previous transfor- 
mation; i.e. at the Ith measurement the basis may be 
transformed by any U^^'^ acting on the subspace (D^~'. 

Proof. Consider measuring in the basis {\al )} on the 
state resulting from the previous measurement, namely, 

|ar^>i|'5if-"i'^(a^'^;ar^)>2,-,7V. Repeating Theorem 
m on the smaller singlet, one can see that the final 

i i • I I (2)\ 10(^-2)^ (21 (2) (2)n\ 

State IS |a, ')2|5^]v-2 ))3,--,Af- 

In general, if M measurements are 
taken, we have (upto a global phase) 

\<Xni)l- ■ ■ \anM )m\On-M [(^nM,^ '))m+1,---,N 

where n = (rii,-- - ,nAf) is a vector, the elements of 
which are the indices of the vectors of the basis Bj^i 
excluded from the [N — M)-singlet. □ 

It must be noted that in the proof — see equation 
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\s 



{N) 
N 



([T]) — we make use of the property U^'^\S[ 
This holds only when the operates on the space 

occupied by the singlet 15^^); if one were to put this 
state within a larger space, a separable unitary opera- 
tion on the whole n-partite space would not (in general) 
give this invariance. Thus in order to iterate the proof, we 
make the restriction that at each successive measurement 



the dimension of the unitary transformation decreases by 
one. So, at the Ith measurement, one would be able to 
transform the basis by a unitary U'-^\ such that there is 
a submatrix operating on TV — Z levels (which levels are 
operated on depends on the previous outcomes), and '1' 
on all diagonal elements corresponding to the remaining 
levels. In fact, this restriction may be lifted slightly. It is 
well known [llj that any dxd unitary matrix can be writ- 
ten as a product of two- level unitaries: Ud — Vi - ■ - Vk, 
where k < d{d — l)/2. This then gives 



Ud(g>Ud^iVi(g>Vi)---{Vk(g>Vk), 



(2) 



allowing one to use at measurement the subset oi d x d 
unitaries that factorise such that the two-level matrices 
{Vi} either act within the subspace supporting the sin- 
glet, or its complement {i.e. the matrix element linking 
the singlet subspace to the rest of the space is zero) . This 
is because those factors operating on the complement act 
as the identity on the singlet, and those operating on the 
singlet subspace have the rotational invariance property. 

But what would happen if one measured in an arbi- 
trary basis? To answer this, consider the effect of oper- 
ating on a 2-singlet with a more general unitary operator, 
written in the form Those factors Vi (S) Vi that do 
not operate on any part of the singlet subspace can be 
removed successively from the right of the operator ([2]) , 
until the right-most factor links one of the levels of the 
singlet with a different level. Suppose the ith factor in 
^ operates on levels j,k in a basis {|0}f=ii a-nd the 
whole operator Ud Ud is applied to a singlet of lev- 
els j,m. Then (Vi » Vi) • • • {V <^ V){\jm) - \mj)) = 
{Vi ® Vi) • • • iV-i ® V-i)iiV\j))\m,) - \m)iV\j))). This 
remains a singlet, albeit in a different subspace. In gen- 
eral, when an n-singlet is operated on by a separable 
unitary Uf"' (where d > n), the state remains an n- 
singlet, but (in general) lies within a different subspace 
of the ^"-dimensional n-partite Hilbcrt space acted on by 
uf". Thus we conclude that for any measurement, we 
still obtain a singlet, but for those measurement bases not 
satisfying the constraints given above, the singlet moves 
into a different subspace. 

Entanglement properties: Though there is a plethora of 
entanglement measures for many-body systems, perhaps 
the most physically-relevant is the amount of entangle- 
ment one can establish between spatially-separated sub- 
systems for subsequent use in networked quantum com- 
putation. Indeed, it was with this motivation that local- 
isable entanglement was introduced 0] and a maximally- 
entangled state of a 2 X 2 system has been shown to 
be establishable this way. However, any realistic dis- 
tributed quantum computation scheme would require a 
larger amount of more entanglement. Using the above 
statements, we can now consider the role of A'^-singlets 
in entanglement distribution. Suppose that Alice and 
Bob each hold n qudits of an A^-singlet. By measuring 
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FIG. 1: Two subsystem localisable entanglement. Here, Al- 
ice and Bob have access to the boxed qudits. By performing 
arbitrary measurements on the other qudits (in random di- 
rections), they can establish a 6-singlet between them. 

the remaining qudits, they can establish a 2n-singlet be- 
tween them. This is equivalent to sharing a maximally- 
entangled state of dimension D x D (where D :— {^^)) 
or log2 D EPR singlets, since the 2n-singlet can be writ- 
ten in a Schmidt decomposition of (^^) terms of equal 
amplitude (see Fig. [T|), and upto local unitaries this is 
equivalent to J2iLi N)!*)/'^/^- Thus the localisable en- 
tanglement in our model (when the notion is generalised 
to two subsystems) can be much larger in comparison 
to the models studied so far. Furthermore, by taking a 
Schmidt decomposition across an arbitrary partitioning 
of qudits, one can see that the block entropy of this state 
is log2 (^) . This is in contrast to the usual behaviour 
for gapped systems, where the entropy is usually propor- 
tional the block's 'area', and to the case of gapless spin 
chains in ID where it is proportional to log2 L. We fi- 
nally state without proof that the A^-singlet is genuinely 
iV-partite entangled (i.e. the state is non-separable with 
respect to all possible partitionings of qudits, and it is not 
producible by fc-partite entanglement, for all k < N). 

Persistency of entanglement: An important considera- 
tion from a practical point of view is the ease of destroy- 
ing the entanglement present. This can be quantified by 
the persistency of entanglement, the minimum number of 
local von Neumann measurements that will disentangle 
the state [10]. This may be used to model an environ- 
ment that interacts with a system through random, local 
measurements. Since the establishment of successively 
smaller singlets is independent of both the basis choice 
and outcome, it is easy to see that — 1 measurements 
always need to be performed to completely disentangle 
the state (there is no way to optimise the basis to reduce 
this number). Thus these states have the highest possi- 
ble persistency of an A^-partite state (for comparison, the 
multi-qubit cluster, GHZ and W states have persistencies 
[N/2\ , 1 and - 1, respectively). 

Physical preparation of states in ID: The permutation 
Hamiltonian may be obtained as an effective Hamiltonian 
in a certain regime of the d-species Hubbard model: 

d 

i,j (T—1 i a^(j' 

where i, j are summed over some network, V^^a' is the 
interaction between species a and cr', t^- the hopping in- 
tegral between lattice sites i and j for species cr, and 



c^j is the creation operator for a at lattice site i (ei- 
ther fermionic or bosonic). Fermionic: In the regime 
of strong, isotropic couplings [U 3> iy ), with l/d-fUhng 
regime, this (following a similar perturbation calculation 
to Fradkin and Emery (T^ for the two-level case) 
becomes H' = Y.'L,f3=iJ2i.j Jijcl^cpiC^^Caj, exactly the 
permutation Hamiltonian required with Jij = Atf^/U 
[26|. Bosonic: the permutation Hamiltonian is realis- 
able in the l/d- filling regime with careful tuning of the 
species-dependent tunneling integrals and the addition of 
a hard-core interaction 

This is one of the most widely-studied condensed mat- 
ter models, and has been shown to be realisable with 
ultracold atoms in optical lattices. Two- level atoms at 
1/2-filling realise the Heisenberg model [l3l, and replac- 
ing these with atoms with d degenerate levels (e.g. d 
hyperfine levels 15] where d = 10 is routine with 
[ig ] and much higher d is possible with Er [l7[) realises 
the permutation Hamiltonian. To prepare the singlet, a 
d-site optical lattice is first loaded with d-level atoms in 
the translationally-invariant state jl)®'* (with all tij 0). 
An inhomogenous B-field coupled to the hyperfine levels 
is applied to select this as the unique ground state within 
the 1 /d-fiUing regime and to allow parity to be broken in 
the subsequent evolution. The ratios tij/U are then in- 
creased adiabatically (the standard method for ground 
state preparation in optical lattices [3]) to reach finite, 
but small tij <C U (still in the Mott insulator regime); 
similarly the _B-field is slowly tuned to zero. The system's 
state evolves to the final ground state: the singlet. Other 
potential phwical implementations include spin ladders 
and tubes [l9|, arrays of quantum dots with both spin 
and orbital levels [6|. To date, little progress has been 
made in realisation, although entangled photons [1^ and 
cavity QED [2l|, [23] have been considered. 
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